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Abstract. We study the behavior of massless Dirac particles in the vacuum C- 
metric spacetime, representing the nonlinear superposition of the Schwarzschild 
black hole solution and the Rindler flat spacetime associated with uniformly 
accelerated observers. Under certain conditions, the C-metric can be considered 
as a unique laboratory to test the coupling between intrinsic properties of 
particles and fields with the background acceleration in the full (exact) strong-field 
regime. The Dirac equation is separable by using, e.g., a spherical-like coordinate 
system, reducing the problem to one-dimensional radial and angular parts. Both 
radial and angular equations can be solved exactly in terms of general Heun 
functions. We also provide perturbative solutions to first-order in a suitably 
defined acceleration parameter, and compute the acceleration-induced corrections 
to the particle absorption rate as well as to the angle-averaged cross section of the 
associated scattering problem in the low-frequency limit. Furthermore, we show 
that the angular eigenvalue problem can be put in one-to-one correspondence 
with the analogous problem for a Kerr spacetime, by identifying a map between 
these “acceleration” harmonics and Kerr spheroidal harmonics. Finally, in this 
respect we discuss the nature of the coupling between intrinsic spin and spacetime 
acceleration in comparison with the well known Kerr spin-rotation coupling. 


PACS number: 04.20.Cv 


1. Introduction 

The nature of the interaction between intrinsic spin and acceleration has been 
extensively investigated in the literature from both classical and quantum perspective, 
in view of possible violations of the equivalence principle between inertial and 
gravitational masses as well as of parity and time-reversal invariance in the 
gravitational interaction (see, e.g., Ref. [1] for a review). Several studies of the Dirac 
equation in uniformly accelerated reference frames and gravitational fields have shown 
that a direct coupling of spin with linear acceleration does not arise [H ini H] . This is 
in agreement with the lack of observational evidence in favor of such a coupling. In 
contrast, there exists a coupling of intrinsic spin with rotation, which is independent 
of the inertial mass of the particle [S] . For instance, in the case of photons the helicity- 
rotation coupling is responsible for the phenomenon of phase wrap-up, which has been 
tested with high accuracy via rotating GPS receivers [5]. Furthermore, in the case of 
neutrinos this kind of interaction may lead to a helicity flip [7] . 

In this paper we study the behavior of massless spin-^ particles in the vacuum 
C-metric background, which represents the exterior gravitational field of a uniformly 
accelerating Schwarzschild black hole, under certain conditions. We refer to Ref. [5] 
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for an exhaustive review of its main geometrical and physical properties as well as 
its historical background. By using spherical-like coordinates, the Dirac equation 
is separable into one-dimensional radial and angular parts. The radial and angular 
equations are coupled via the angular eigenvalue. We consider first the angular part 
as an eigenvalue problem, discussed here to leading order in perturbation theory, 
following the approach of Press and Teukolsky [9] . We provide the corrections to the 
angular functions as well as to the associated eigenvalues to first order in a suitably 
defined acceleration parameter. The radial equation is then solved perturbatively 
up to a certain post-Newtonian order. The particle absorption rate as well as the 
angle-averaged cross section of the associated scattering problem are also computed 
in the low-frequency limit. Noticeably, both radial and angular equations can be 
solved exactly in terms of general Heun functions, which are but of poor practical use 
and only of formal utility. Finally, we discuss some features of the coupling between 
helicity and acceleration. 

Units are chosen so that G = 1 = c = fi and the metric signature is —2. 


2. Dirac equation in the C-metric 


The C-metric belongs to a class of degenerate metrics discovered by Levi-Civita [10] 
and takes his name from the classification of Ehlers and Kundt m- It represents 
the exterior field of a uniformly accelerated spherical gravitational source and it 
can also be thought of as a nonlinear superposition of the (flat) Rindler spacetime 
associated with a uniformly accelerated family of observers and the Schwarzschild 
solution for a static black hole [nm. The line element has been expressed using many 
different coordinate systems, those more suitable to study their geometrical properties 
[iKiniiiiiiTKiaiiij. We will adopt Schwarzschild-like coordinates x°‘ = {t,r,0,(j)), 
in terms of which the line element writes as 1191 


ds^ = 


1 


Q{r) 


dr - 


dr^ - 


■ de'^ -r'^p{e) sm"^ ed(j)^ 


( 2 . 1 ) 


U2(r,6») V Q{r) P{0) 

where Q{r) = r{r — 2M)(1 — P{9) = 1 — 2MAcos9 and D.{r,9) = 1 — Ar cos 9. 

The constants M > 0 and A > 0 denote the mass and acceleration of the source, 
respectively. The metric (ED) then reduces to that of a Schwarzschild black hole for 
^ = 0 and to the flat spacetime in uniformly accelerating coordinates for M = 0. 
In this coordinate system, r is constrained to the interval 2M < r < 1/A between 
Schwarzschild and Rindler horizons. Furthermore, one should require P{9) > 0 to 
preserve metric signature, i.e., cos0 < 1/{2MA), thus excluding an entire conical 
region around the positive z—axis. However, if one limits the acceleration parameter 
so that MA < 1/2, then P{9) > 0 for all 9 G [0,7r]. Hereafter we will assume 
this condition to hold. Unavoidably, a conical singularity still exists, since P(0) = 
1 — 2MA ^ 1. The latter can be removed by limiting the range of allowed values of 4> 
in the interval (j) = [0, 27r/P(0)]. No further constraint comes from the conformal factor 
U, which is everywhere positive between the Schwarzschild and Rindler horizons. 

Let us introduce the following Kinnersley-like null frame [20] 


2^2 


I = 


m = 




Q 

n 

VTr 


a 

Ot + -^o, 

[y/Pde ■ 


1 


n = 


a ^ a 

dt - 7^0, 


\fP sin 9 


d. 


( 2 . 2 ) 
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By using standard notations and conventions of the NP formalism (see e.g., Ref. 
m) we analyze below the dynamics of Dirac particles (with special attention to the 
massless case). The wave function associated with such particles is written in terms 
of a pair of spinors and (spinor indices are denoted by capital letters and run 
from 0 to 1), whose components are often indicated as = Pi, = F 2 , = Gi 

and Q° = —G 2 , resulting in the following equations 


{D + e — p)Fi + (d* + TT — a)F 2 = ip*Gi, 

(A + /r - 7)^2 + ((5 + ,5 - t)Fi = 

[D + e* - p*)G 2 -{5 + M - a*)Gx = ip,F 2 , 

(A + - 7 *)Gi - {5* + r - t*)G 2 = ip.Fi ; (2.3) 

here /r, is the particle mass and D = A = n^df^, S = m^i9^ and <5* = 

denote frame derivatives. The nonvanishing spin coefficients associated with the frame 
(lO) are given by 


P = 


D (D — r Or) 


r = —TT = —llg 


Q p 


P = 


O sin 6 Pg + 2 P cos 6 


V2r ’ ^ 2r2 O^ 

1 rOQr - 2 rQOr - 2 Q O 


a=^ne-l3, 


\/2r 


r/Psin6» ' ^ 4 r^D 

where Xr = drX and Xg = dgX. The only nonvanishing Weyl scalar is ip 2 = 
-Mn^/r^. 

The existence of a time-like Killing vector dt and a rotational one in the C- 
metric allows us to assume a wave function with the customary dependence , 

It is useful to define the following radial and angular differential operators as 


(2.4) 


d iujr'^ Qr 

"" ar <3 + " 0 ’ 


I>!. = (I>.r = S:+ <3 


Qr 


(2.5) 


and 


„ d mcscO ^ ^ , d to esc 0 „ , , 

Cn = -^4 - P - ^"n-cote, Cl^Cni-m) = — - - -hncot6», (2.6) 

so that D = O^Do, A = — 6 = and 6* = The following 

rescaling of the quantities Pi _2 and Gi _2 

„pl/4 

h = -^F,, 




9i = —- Gi, and 32 = G 2 , (2.7) 


O 02 02 

allows for considerable simplifications in Eqs. (1^ . which become 

P 


O 


VQVofiF 


■^ 1/2 = 


ip*rgi 
0 ^ 


^vlh-VWc\h= -2'j^ 

2 ll 

\ Z 2 \l 

^/QVlg, + V^C.g2 = - 2 ^^^ 


( 2 . 8 ) 





















Massless Dirac particles in the vacuum C-metric 


4 


Let us consider the case of massless particles, i.e., /i, = 0. Eqs. (j2.8l) can then be 
solved by separation of variables. In fact, defining 


/i(r,0) = / 2 (r, 0 ) = i?+i(r)5+i(0), 

gi(r,0) = E_^i(r)S_i(0), g 2 (r,d) = -^E_i(r)S+i(d), 


(2.9) 


we get the radial equations 

\/QDoR_i = XRx, y^vlRx = XR_i, (2.10) 

and the angular equations 

VpCiSi =-XS_i, Vpc\S_i=XSi, ( 2 . 11 ) 

22 2 ■222 

where A is the separation constant. Note that the complex conjugate angular functions 
51 1 satisfy the same equations (j2.11l) as 5+i. 

It 2 2 


3. Angular equation 


By applying the differential operator C\ to the first equation of Eqs. (12.111) and 

2 

using the second one to eliminate 5_ 1 , we get the following second-order differential 
equation for the angular function Si 

d^^i / P, \ d5i 


0 = 


d6'2 

1 

“ 2 


g+cot0 


d6» 


1 - 


2 A 2 


Pe 


-{2m — PcosB) 


(2m -b P cos 0)2 


Si 


(3.1) 


P ' 2P2sin0'""” " ’ 2P2sin2 0 

The corresponding equation for S_i can be obtained from Eq. simply by 

replacing m —>■ —m. Both cases s = can be handled together by introducing 
the following equation 


1 


sin 6 d0 
where ^5(0) = P^/'^S^x and 
l + E-s"^ 


sin0- 




d0 J 


+ ^(ang)(0)«5(0) = 0 , 


(3.2) 


^(a„g)( 0 ) = 


1 

p2 


[s -b (m — 2sMA) cos 0]2 


sm 


-b (m -b sMAY -b (1 — s^)(l — MAcosOY — M‘^A‘ 


(3.3) 


where the separation constant has been replaced by A^ = P -b s^. We choose the 
following normalization for the angular functions 

y'd0^|,5(0)p = l, (3.4) 

also implying 

y'd0^|5±i(0)|2 = l. (3.5) 

Eq. (13.2|) represents a Sturm-Liouville eigenvalue problem for the separation 
constant E. We will solve it perturbatively in the next section to first-order in the 
acceleration parameter ry = 2MA, following the approach of Press and Teukolsky [^. 
We will then provide an exact solution in terms of Heun functions. 
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Following Press and Teukolsky [3] , Eq. (13.21) can be written as an eigenvalue equation 
involving the sum of two operators, i.e., 

{no + vni)sS = -EsS, (3.6) 

to first-order in the acceleration parameter rj, where 


Ho= - 


1 d 


sin6ld6» 




-I- -I- 2ms cos 0 


sin^ 9 


Til = ms + L cos 6 — 2m 


s + m cos 9 


sm 


(3.7) 


with L = l{l 1). Let 

E = L f]Ei , s^lm — s^lm 'n s^It: 


(3.8) 


where sYim{9) denote the spin-weighted spherical harmonics (SWSH), with I = 
|s|, |s| -I- 1,... and —l<m<l. We then have 

Hq sYlm — L sYlmj 

cL{L - L')sYi,m =-{%!+ E^),Yl^ . (3.9) 

Multiplying (13.91) 9 by sin^ ^ sY{i^ and integrating over the solid angle dO = sinddddc/i 


gives 


C\^[L — L')( sin^ = — {El tos)( sin^ 9)^ — L( cosd sin^ 9)^ 


(3.10) 


-I- 2ms6ii> + 2to^( cos6()j , 
where the following notation has been used 

= {sl'm\X\slm) = J dn sY;^X sYi^. (3.11) 

can be computed by using the relations between the SWSH, the Wigner 
rotation matrices and the Clebsch-Gordan coefficients (cf. [9] and references therein). 
Multiplying Eq. (011)2 by sin^ integrating over the solid angle gives instead 

Cl^{L — L')( sin^ j = — {El + ms)( sin^ 9)^ — L( cosdsin^ 0)^ 

+ 2ms + 2m^( cos 9)^. 

If Z = /', one gets the following solution for Ei (from Eq. (13.121) 1 

El = —ms -I- — \ . \—L( cosdsin^ 9), + 2ms + 2m}{ cosd),!. 

(sin^d)^ ' ^ 

If Z Z', one gets the following solution for (from Eq. (13.101) 1 


cL = 


L-U 


-{El -I- ms) 


(cosdsin^d) A cos 9) 


t^-L 


2m^ 


{smH}^, (sin^d);, (sin"^);, 

Therefore, we need to calculate the following quantities 
(sZ'm| cos0|sZm) = (ZI ?7i0|Z'm)(ZI — s 0|Z'— s), 

(sZ'm| cos^ 0|sZm) = ^Z'+i (Z2 mO|Z'm)(Z2 — sO|Z' — s), 

(sZ'm| cos^ 0|sZm) = {llmQ\l'm)(ll — s0\l'— s) 

-|-| (Z3mO|Z'm)(Z3 — sO|Z' — s)) , 


(3.12) 

(3.13) 

.(3.14) 


(3.15) 
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where (jij 2 Wim 2 | are the Clebsch-Gordan coefficients, which are nonzero only 
if M = TOi + m 2 , with J = \ji — j 2 \T--,ji + j 2 and M = — J, The 

explicit expressions for the above quantities are shown in [Appendix A[ The first-order 
corrections to the energy eigenvalues turn out to be (see Eq. (IA.6I) 1 

2[P{1 + 1)2 + 5m2] - l{l + 1)(7 + 2m2) 


El = ms- 


[i-m 


(3.16) 


- 2)[1(Z-I-1)-I-m2] 

are given by Eq. p.lSp . The first-order 
corrections to the eigenfunctions follow directly from Eq. (13.81) . For example, for 
/ = 1/2, 3/2, 5/2 they result in 


whereas those to the coefficients 


+ CL.nr 


, 5 . = CaL.n,. + Ci^sY. + 


Cl sYi^ 


— CIsYx + 


Cl sY3 

2 


Cl sYv 
Sm® 2 


+ cL.n, 


■Cl La 

2 


+ CLLii 


(3.17) 


where the coefficients are listed in Table [TJ The behavior of the eigenfunctions 
Li = is shown in Figs. [T] and [2| for a given value of the acceleration 

parameter for the lowest modes Z = 1/2 and I = 3/2, respectively. 


Table 1. The numerical values of the coefficients needed to compute the 
first-order corrections II3.17II to the angular eigenfunctions are listed below. 



^Im 

1 

m 

V - ^ 

1' - ^ 

^ 2 

V - ^ 

^ 2 

V - ^ 

^ 2 

V - ^ 

^ 2 


1 


\/2 

23\/3. 

9 



1 

2 


80 

1080 

800 



2 

1 


V2 

23V3s 

9 




2 


80 

1080 

800 




3 

0 


39 

VlOs 

11 V 5 



2 


154 

1680 

882 



1 

53 V 2 


71\/6 

39V2s 

33vT0 


3 

2 

280 


1512 

8960 

2450 


2 

1 

53 V 2 


71 Ve 

39\/2s 

33vTo 



2 

280 


1512 

8960 

2450 



3 

0 


39 

VTos 

llV^ 



2 


154 

1680 

882 



5 

0 

0 


145\/6 

II5V2IS 

65\/l4 


2 


792 

31248 

6804 


3 

0 

43 


113%/l0 

bVbs 

65Vi35 


2 

462 


3960 

3024 

6156 


1 

169v^s 

467^6 


323 V3 

17vd5s 

325 V 2 

5 

2 

1080 

3024 


2160 

7560 

5832 

2 

1 

169v/3s 

467\/6 


323 V 3 

nVTEs 

325\/2 


2 

1080 

3024 


2160 

7560 

5832 


3 

0 

43 


113vT0 

5-\/5s 

65\/35 


2 

462 


3960 

3024 

6156 


5 

0 

0 


145\/6 

115v^. 

65\/T4 


2 


_Z22_ 

31248 

_eS24_ 


3.2. Exact solution 


Eq. (13.2|) can also be solved exactly in terms of Heun functions. In fact, suitably 
rescaling the function sS as 

2 mr] mrj 


5 = yPL'=+(z-l) 


k- 


2h — s -(-!-(- 


rf - 1 


2k± = s — 


?7 ± 1 


(3.18) 
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(a) (b) 

Figure 1. The behavior of the angular eigenfunctions S,i is shown for 1 = 1/2 

^ 2 

and a value 7] = 0.3 of the acceleration parameter, as a function of the polar angle 
9. Panels (a) and (b) correspond to m = —1/2 and m = 1/2, respectively. Dashed 
curves are the corresponding (unperturbed) SWSH (77 = 0). 


in terms of the new variable z = cos^ Eq. (03) becomes 

^^ 1 / (7 S e \dy ap-q ^ 

dz'^ \z z — 1 z — a) dz z{z — 1 )( 2 ; — a)^ 

which is a General Heun equation |23j in standard form with solution 


y = ciHeunG[a, q, a, /3, 7 , <5; z] 

+ C 2 Z^“'''HeunG[a, q + (1 — 'y){a5 + e), /3 + 1 — 7 , a + 1 — 7 , 2 — 7 , d; z] , 


(3.19) 


(3.20) 


and parameters 

2r]a = 1 + ry, 2r]q = (1 + 2r])s^ + (1 + 37 )s + q — E , 

a =l + s, /3=l + 2s, 7 = 1 + 2k+, <5=1 + 2fc_, (3.21) 

with 7 + (5 + e = a + /3 + l and a ^ 0,1. This equation has four regular singular points 
at Zi = {0,1, a, 00 }. In fact, the coefficient of y' in Eq. (13.191) has a singular part 
diverging as Ai[z — Zi)~^, whereas that of y behaves as Bi{z — Zi)~'^ + Ci{z — Zi)~^, 
with at least one of the coefficients Ai, Bi and Ci nonvanishing. Furthermore, the 
General Heun function is such that 

HeunGla, q, a, j3, 7 , d; 0] = 1, —HeunGla, q, a, /3, 7 , d; 0] = — , 

dz 7 a 

and its expansion around z = 0 starts with 

HeunGla, q, a, f3, 7 , d; z] = 1 + —z 

7 a 

^ q + qaS + q'/a - q6 + qa + qjd + q^ - 2 , 3 ^ 

7(7 + l)a 2 z+UGj. 


(3.22) 


(3.23) 
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(a) (b) 




(c) (d) 

Figure 2. The same as in Fig. [T] but for the I = 3/2 mode. Panels (a) to (d) 
correspond to m = —3/2, m = —1/2, m = 1/2 and m = 3/2, respectively. 


4. The radial equation 


Let us consider the radial part of the Dirac equation. Combining Eqs. (|2.10l) . we get 
the following second order equations for the radial function R_i 


which reads 
d^R 1 


iQV\Vo - X^)R_i =0, 
2 2 

Qr d-R-i 


dr^ 2Q dr Q 


A" - 




+ 2iujr 1 — 


rQ'r 


R_i = 0, 


(4.1) 


(4.2) 
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whereas i?i satisfies the complex conjugate equation. Both cases s = ±i can be 
handled together by introducing the following equation [52] 

= 0 , (4.3) 

with 

2 4 

^(rad)(^) = — 2 ryl^(r — M)(l + s)(l + 2 s) H— 

- (dw - idiv) - + ”<1 + 2.). (4.4) 

and the new radial function is such that _iTZ = R_i and \/Q iTZ = Ri. 


4 . 1 . Perturbative solution 


Let us consider first perturbative solutions to the radial equation (|4.3I) for small values 
of the acceleration parameter 77 = 2MA. Since A enters this equation quadratically, 
to first-order in ry the radial equation maintains the same form as in the Schwarzschild 
case (77 = 0), differently from the angular equation. Nevertheless, a dependence on the 
acceleration parameter still remains through the separation constant, which is such 
that = E + s'^, and the energy eigenvalues E = L + pEi are affected to first-order 
in 77 according to Eq. (IA. 6 I) . 

We list below the first few terms of the corresponding PN solution (M ^ e^, 
A ^ and uj ^ e, where e = 1/c is the inverse of the speed of light) to first-order in 
77 . The two independent solutions are 


s^(in) — s^(in) + VEl s^[in) ) s^(up) — s'^(up) + s^(up) : (4-5) 

where the subscripts “in” and “up” refer to the regularity for small and large values 
of r, respectively. We find 


(0) 


isujr 

7TT 


e — 


{l-s)— + 


isu) 


(I - l)[2l - I - + 


r 4(; + l)J 4(; + l)(Z + 2) 
m2 


4 4 
uj r 


8(/ + 1) 2 ) _ 

'JD (0) _ 'T? I 


r- {I + 1)2 

d + 0(edl , 


Moj r 


and 

— 1^21 + 1)2 
isLor 


- 1 


(/ + 1)2 

21 - 1 


{21 + l)ln 

(Z + l)( 2 Z + l)lnf 


r 
\'^o 


-31-2 


4(1 - s){l - s - 1) 


M 


[{21 + 1){2P - I - 6ls + l + s)\n 


ro 


- {21 + 5)s + 21^ + I - 2] 


4.{l + l)‘ 


:{2l - 1 - 4s)[(/ -b 1)(2; -b 1) In — - 3Z - 2] 


ro 


(4.6) 

























Massless Dirac particles in the vacuum C-metric 


10 


isuj^r^{2l — 1)(2Z — 1 — 4s) 
16(Z + 1)2(Z + 2)2(Z-s)(Z-s-1) 
-5P- llZ-5] 

21-1 


(Z + l)G + 2 )( 2 Z + l)ln - 


16(Z — s){l — s — 1) [ r2 
+ (-16Z + 24 - 32/2)s - 5; + 12?^ + 12 - 


\/f ^ ( r 

[(/ - 1)(2; + 1)(4/2 - 16/s - 4/ + 8s + 5) In ; ^ 


Muj^ 


[{I + 1){21 + 1)(2/3 - 13/2 - 6/2s -71 + 25/s + 28s + 7) In 


(/ + 1 ) 

+ lO/^s - 87Ps + 50/2 _ 59^ _ 14 _ 2/^ + 43^3] 

+ 8{i + mi + ^y^^^ ~ ^ ^ ^ ^ 

-5/2- 11/-5]] e'^ + 0{e^)) , 

T? (1) _ T? (1) I 

where ro is an arbitrary scale factor. The constant Wronskian 


(4 


W = Q 


s+l 


7^u;Zm(in)^ sT^wZm (up) s'/^cuZm (up) ^ sT^cuZm (in) 


dr ‘ 


(4.8) 


is given by 


W= - 21-1 + isMuj 


1 - r]Ei 


[2l{l + 1) - 7s - 4](2/ - 1)(2/ + 3)(2/ + 1)^ 
16/(/ + 1 + s)(/ — s)(l + 1)(/ — s — 1)(/ + 2 + s) 
4s/(/ + 1) — 7s — 4 1 3 5 2^ 

e^ + 0(e^?72). 


(2/ + l)2/(/ + l)[2/(/ + 1) - 7s - 4]_ 


(4.9) 


4-2. Exact solution 


As in the angular case, the radial equation (14.31) can be solved exactly in terms of 
General Heun functions. In fact, suitably rescaling the function gTZ as 

= wQ-^r - 2Mf^ (1 + Arf^ {Ar - 1)'=^ , (4.10) 


with 


ki = 


2iujM 


k 2 = s + 


iujM 
77(77 + 1 ) ’ 


in terms of the new variable 

r 1 - 2AM 
^ “ W 1-Ar 


fca — — 1 + 


iu>M 
77(77 - 1) ’ 


(4.11) 


(4.12) 


interpolating between z = 1 {r = 2M) and 0 —)> 00 (r = 1/A), Eq. (14.31) can be cast 
in the form of a General Heun equation (13.191) with solution w = w{z) given by Eq. 
(13.201) with parameters 

2 r]a = 77 — 1 , 2 r]q = 77(1 — s 2 ) + E, 

a =l + s, /3 = —(1 + 2 /C 3 ), y = l —s. 


S = l-s + 2ki. (4.13) 
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In order to study the asymptotic behavior of the radial functions at the Schwarzschild 
and Rindler horizons it is useful to introduce the tortoise-like coordinate r, defined 
by dr^f/dr = r"^ jQ, i.e., 


r. 


—--—- In Vl + Ar -^—- 

A{1 + 2AM) A{1 - 2AM) 

2M , z' r \ 

+ TV4M535‘"(2M-d+““‘- 


In \/l — Ar 


(4.14) 


It is a single-valued function of r, because dr/dr^ is always positive, and is such that 
r, —>■ —oo when r —>■ 2M, and r* —>■ -|-oo when r —>■ XjA. 

By introducing the scaling sAR, = r~^Q~'^H, in terms of the new variable r* the 
radial equation (j4.3l) can be transformed into the one-dimensional Schrodinger-like 
equation 


+ VH{r) = 0 , 


with potential 


V = 


Q 




(rad) 


2 Q Qr 


--,IQ 


+ 


s Q 


2 \ 1 


The asymptotic form of the radial equation as r - 
d^ 2 

-^H{r) -I- [w -I- isKA^ H{r) = 0 . 
drf 


2 Q 

\/A(r If — I 


oo) is 


(4.15) 


(4.16) 


(4.17) 


where ka = ^(1 — 2MA) is the value of the surface gravity at the Rindler horizon. 
For small values of A we then find 

[oj + isKA]^ = + 2iAsuj + 0(A'^) . (4-18) 


The solution is then given by 


H 


^±i[bJ + i8KA\r» 


QT^/2g±i<^r. 


(4.19) 


implying that AR ~ and AR ~ for outgoing and ingoing waves, 

respectively. 

On the other hand close to the Schwarzschild horizon r —>■ 2M (r* —>■ —oo), the 
asymptotic form of the radial equation becomes 

d^ 2 

—+ [a; - H{r) = 0, (4.20) 

where k+ = {1 — AM‘^A^)/(AM) is the value of the surface gravity at the Schwarzschild 
horizon. For small values of A we then find 

[oj — isKj^f' = (oj — 'isK+)^ + O(A^). 

Therefore, the asymptotic solution is given by 


(4.21) 

(4.22) 


implying that sR (outgoing waves) and sR (ingoing waves). 

The radial equation (14.311 is associated with a one-dimensional scattering problem, 
once suitable asymptotic boundary conditions are imposed (see Ref. [5]). An incident 
wave traveling towards the Schwarzschild horizon with a given amplitude will be 
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partially reflected by the potential barrier (reaching then the Rindler horizon), and 
partially transmitted across the black hole horizon. This situation is realized by the 
conditions 

T —>■ 2M 


.JZ 






i?e*' 


1 

A 


(4.23) 


where i? and T are reflection and transmission coefficients, respectively. 

5. Dirac current for massless particles 

Let us study now the Dirac current for massless particles, i.e., 

J>^ = V2 a%, {P^P^' + Q^Q^') , 
where cr^^, are the generalized Pauli matrices defined as 


^AB' — 


V 2 




(5.1) 

(5.2) 


In terms of the rescaled functions ft and gi it reads 

j = [d/ip + \ 92 ?)^ + (I/2P + \gi?)n + (/1/2* - g2gl)m 

+(/2/r - gigD'fh] ■ (5.3) 

and is conserved, i.e., = 0. One can then define the time rate of the number of 


particles entering a r = const, hypersurface as 

^ = -1 


(5.4) 


where g is the determinant of the spacetime metric, so that ^/—g = sin 0/0'^. 
Substituting the assumption (12.91) into Eq. (15.3|) we find 


J<l> = -2JP- 


(5.5) 




r'^y^QP sin 9 

The angular equations (|2.1ip imply 

sin 0 Im(S'_iS'?) = const., (5.6) 

which we set to zero. Similarly, the radial equations (j2.10p yield 

|i?|p-|i?_iP = C, (5.7) 

where C is an integration constant. Therefore, the Dirac current turn out to be given 

by 


J = 


2 QVp 


(l^i p + |5_i|2) (|i?, |2; + |i?_. pn) = J+ + J_ , (5.8) 


so that J^Jfi = 2J+ • J_ > 0, whence J is in general a time like vector. Since the 
only non-trivial spatial component of is J*", we analyze the particle flux through 
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spherical surfaces. If we assume the normalization (13.51) of the angular functions, then 
the conserved net current of particles becomes 


1 dN 
27r dt 


= C. 


(5.9) 


Turning then to the scattering problem discussed in the previous section, one can 
calculate the rate at which particles enter the Schwarzschild horizon per unit time, 
i.e., 


2 tt \ dt ) r_^2M ^ ^ 2 /, 


>2M 

The asymptotic behavior (14.231) of sTi- for r —>■ 2M implies 

so that 


(5.10) 

(5.11) 


which is always 
The net current 


^(^) =m= 


dN 


(5.12) 


>2M 


positive. Therefore, superradiance cannot occur in this case [22 
of particles crossing the Rindler horizon is instead given by 


1 

27r V dt 


1-\R\\ 


(5.13) 


Let us evaluate the absorption rate of particles from the Schwarzschild horizon by 

using the perturbative solution of the previous section. For vanishing acceleration the 

transmission coefficient per unit amplitude of the incident wave must reduce to that 

computed by Page in the case of a Schwarzschild black hole, i.e., FT^'" = 

2 

to the leading order approximation (i.e., for Muj 1, which is the range of validity 
of our approximate solution, and I = 1/2). In order to calculate the transmission 
coefficient |T| = |i?i | = \^/Ql'R\ we need the value of the radial function at r = 2M. 
The latter will be a superposition of ingoing and upgoing solutions (14.51) . i.e., 


l7^ = iC'(in) l7^(in) + iC'(iip) l7^(up) , (5.14) 

where the coefficients i C'(in) and iC'(up) both depend on I and w. Their explicit form 
can be determined by using standard techniques (see, e.g., Refs. [131|5B]). However, 
to first order in the acceleration parameter the correction turns out to be proportional 
to the Schwarzschild value 

2 


Fi = F? 


schw 


(i - |f;i) = M^oj^ (^1 


:mr] 


(5.15) 


where m = ±1/2. 

Furthermore, one can evaluate the low-frequency (angle-averaged) absorption 
cross section 

cr 1 (uj) = F 1 = 2'kM‘^ , Mui <C 1, 

Im 

so that no corrections to the Schwarzschild result arise to that order. 


(5.16) 
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Since neutrinos possess only one state of polarization, they can be described in terms of 
only two nonvanishing spinor components m- Therefore, limiting our considerations 
to left-handed neutrinos (corresponding to the spinor), the Dirac current (15.3|) 
becomes 


J = + l/2p^ + /i/2"2 + f2fim]. 

Substituting then the ansatz dSll) leads to 


(5.17) 


J* 


r 




2QVP 




Re I R 


(«-*« 4 ) 


S_iSX. 


r'^y/QP sin 9 


Im 


(«-4^i) 


S_iSX. 


(5.18) 


Hence, in contrast to the case with and both nonzero, if only one of the chiral 
components of the spinor field is considered, then J is a null vector with nonvanishing 
angular components. 


6. Comparing acceleration and rotation effects 


We have solved in Section 13.11 the eigenvalue problem associated with the angular 
equation given by Eqs. (EHi-dSZl), to first order in the acceleration parameter ij. 
The first order corrections to the energy eigenvalue are given by Eq. (13.161) . It is 
interesting to compare the present analysis with the corresponding one for massless 
Dirac particles in a Kerr spacetime, in order to make a parallel between these two 
complementary situations of uniform background rotation and acceleration. 

Let us briefly recall the Press and Teukolsky result in the Kerr case to first- 
order in the rotation parameter a. The eigenvalue equation reads as TIS = —ES, 
where 


n = 


1 d 
sin 6 d9 



m'^ s^ 2ms cos 9 
sin^ 9 


2aujs cos 9, 


( 6 . 1 ) 


S denotes the spheroidal harmonics and 


E = L — 2auj 


s'^m 


( 6 . 2 ) 


is the energy eigenvalue. Eq. (lOl) shows that for a given value of auj the energy 
depends quadratically on s, whereas for the C-metric it is linear in spin for fixed rj 
(see Eq. (13.161) 1. so that in the latter case the energy levels turn out to split for 
particles with different spin. 

We will show below that quite interestingly the two eigenvalue problems can be 
cast exactly in the same form, allowing to make the comparison easier. By introducing 
the new spin variable s = s -|- mrj, the Hamiltonian operator (Irfl) can be written as 


n = 


1 d 
sin 9 d9 


sin0 


d9 


-I- -b 2ms cos 9 


sin 


r]{ms -b L cos 9), 


(6.3) 


so that the eigenvalue equation (1X71) becomes H sS = —E gS. Notice that L (or 
equivalently 1) can in turn be suitably replaced by L —^ L-\-r]Li without changing the 
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unperturbed Hamiltonian. The constant term ms in the first order operator (lO) can 
then be reabsorbed by a redefinition of energy, i.e., E ^ E rjms. This leads to a 
new eigenvalue problem H gS = —E gS, with 


n = 


1 d 
sin 9 dO 



m? + + 2ms cos 9 

sin^ 9 


+ rjL cos 9 , 


(6.4) 


which formally reproduces exactly the Kerr problem (IQ) with the identifications 
s s and —2aujs O ijL, giving in turn to the spheroidal wave equation a broader 
meaning. The energy eigenvalue (13.161) can also be cast in the form (16. 2p by using the 
freedom in the choice of Li. As a result both angular eigenvalues and eigenfunctions 
in the C-metric case have a one-to-one correspondence with the Kerr energy spectrum 
and spin-weighted spheroidal harmonics. A similar circumstance was observed in 
Ref. |28j . where the comparative analysis of the Teukolsky master equation in the 
Schwarzschild and Taub-NUT spacetimes revealed the (exact) correspondence of the 
spin-weight with a new spin variable depending on the NUT parameter £ according 
to s O s -I- 2ui£. The “quantization” property for the newly defined spin-weighted 
parameter s (which must be a half-integer) thus implies a “quantization” property 
for the acceleration parameter r]. This genuinely (unexpected) new result for the 
acceleration has no counterpart with the rotation, for which no “quantization” rule 
applies in this context. 

Finally, it is worth to mention that the couplings of the spin-weight with the 
rotation and with the acceleration (both of them inertial effects) have different origin. 
In fact, the spin-rotation coupling (also referred to as Mashhoon effect, see Ref. 
0) appears as a direct coupling (as), whereas the acceleration always enters the 
Hamiltonian through the parameter r] = 2MA, so that any coupling with the spin is 
indirect {AMs), i.e., mediated by the mass of the black hole. This is in agreement 
with the results of Ref. [4] . 


7. Concluding remarks 

We have investigated the behavior of massless Dirac particles in the spacetime of the 
vacuum C-metric, which describes the static region in the neighborhood a uniformly 
accelerating Schwarzschild black hole, under certain conditions. Because of the type 
D character of the solution, the Dirac equation is separable in a suitable coordinate 
system. We have adopted spherical-like coordinates and solved the one-dimensional 
radial and angular equations. We have considered first the angular part. Treating 
the angular equation as an eigenvalue problem, we have computed the first-order 
corrections in a dimensionless parameter associated with the background acceleration 
with respect to the unperturbed Schwarzschild case. The solution for the angular 
eigenfunctions has been obtained as a series expansion in the basis of spin-weighted 
spherical harmonics. The associated eigenvalues turn out to depend linearly in spin, 
differently from the case of a Kerr spacetime to first order in the rotation parameter, 
where the dependence is quadratic in spin. Therefore, the energy spectrum splits 
between particles with different spin due to the spacetime acceleration. Concerning 
then the radial equation, we have explicitly computed the first terms of the post- 
Newtonian expansion to first order in the acceleration parameter. We have analyzed 
the associated scattering problem and computed the corrections due to the acceleration 
to the particle absorption rate as well as to the angle-averaged cross section in the 
low-frequency limit with respect to the corresponding results for the Schwarzschild 
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solution. Furthermore, we have provided the exact solution for both radial and angular 
equations in terms of general Heun functions. Finally, we have discussed the nature 
of the coupling between intrinsic spin and spacetime acceleration in comparison with 
the well known spin-rotation coupling in a Kerr spacetime, by writing the angular 
eigenvalue problem exactly in the same form as in the Kerr case (linearized in the 
rotation parameter). This formal analogy implies a suitable redefinition of both the 
spin-weight s and the quantum number 1. We have thus identified a map relating the 
C-metric angular eigenfunctions and eigenvalues with the Kerr spheroidal harmonics 
and energy spectrum, implying a “quantization” rule for the acceleration parameter. 

Appendix A. Perturbative angular solution: coefficients 

We list below the explicit expressions of the various quantities necessary to calculate 
the first-order corrections to both energy eigenvalues and eigenfunctions given by Eqs. 


and (l-mn . 


The relevant Clebsch-Gordan coefficients to compute ( cos 0) 



(A.l) 


and 



(Zl — s0|/ — s) 


s 


\/K^ + 1 ) 



(A.2) 


leading to 


/ 


'JP — tv? 

21 


V = 1-1 



(cosd);^;, = < 


(A.3) 
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, cos = 


—ms 


^ \/P — 

"""?(/-l)(l + l) ’ 

+ 1) — m? 
2;(Z- 1)(/ + 1) ’ 

\/(r+T)2r^7^ 

Z(Z + l)(Z + 2) ’ 


\/{l + 1)^ — m‘^\/{l + 2)2 — m? 


I' = 1-2 

V = 1-1 

V = I 

I' = 1 + 1 
, I'= 1 + 2 


and 


(cos^O) = < 


4(Z + l)(Z + 2) 
y/P — m^^J{l — 1 )^ — m^^J{l — 2)2 — 


-Sms 


—3ms 


8 l{l-l){l- 2 ) 

\/P — m^yj{l — 1)2 — 
Alii - 2){l - 1){1 + 1) 

y/P — m^{P — 1 — mP) 
8l{l-l){l + l) ’ 

[Z(^ + 1) — 1 — mp] 


2 l{l-l){l + l){l + 2 ) ’ 

, \/(^ + 1)^ — m‘^[l(l + 2) — mp] 


8 l(l + l){l + 2 ) 


Sms 


y/il + l)^ - m^y/il + 2)^ - m'^ 
Al{l + l){l + 2 ){l + 3) 


y/ {I + 1)^ — m"^ y^(Z + 2)^ — \/(^ + 3)^ — m? 


(A.4) 


I' = 1-3 

I' = 1-2 

I ' = 1-1 
l' = l . (A.5) 
Z' = Z + 1 
Z' = Z + 2 

Z' = Z + 3 


8(Z + l)(Z + 2)(Z + 3) 

The first-order corrections to the energy eigenvalues (13.1311 then turn out to be 
2{L^ + bm"^) - L{7 + 2m^) 


El = ms- , ^ ^ , 

{L-2){L + m^) 

whereas those to the coefficients Cj^ (see Eq. (|3.14ll l are 

._3 _ y/P - m^yj[l - 1)2 - m^y/{l - 2)2 - m^{l - 3)(Z -|- 1) 
“ 24(Z- l)2[(Z + 2)(Z + 3) + m2] 

^ 1-2 y/P ~ m'^y/{l — 1)2 — TO^[—2tos(Z -I- 1) -|- £'i(Z — 2)] 

^Im 


(A.6) 


^Im 


CL = 0 , 


4Z(2Z-l)[(Z-2)(Z-l)-hTO2] 
y/P — m^{Z(Z — 1)[(Z -I- 1)^ — 5m^] — 2m^(Z — 5) -|- SmsEi} 
8l{l + l)[l{l - 1) + m"^] 
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;_i_i _ \/(Z + 1 )^ — m?{{l + 1)(Z + 2){P — 5m^) + 2m^{l + 6 ) + SmsEi} 

“ 8l{l + l)[{l + l)il + 2) + m^] ’ 

;_i_2 _ l\/{l + 1 )^ — ni^~^{l + 2)2 — m2[(Z + 3)i?i — 2lms\ 

^ 4(; + i)(2; + 3)[(; + 2)(; + 3) + m2] ’ 

.+3 ^ _ l{l + 4)^{I + 1)2 - mV(^ + 2)2 - mV(^ + 3)^ - 

24(Z + 2)2[(; + 3)(/ + 4)+to2] ■ ^ > 

Finally, the corrections to to the eigenfunctions then immediately follow from their 
definition f|3.8p . 
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